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NOMENCLATURE 


DESCRIPTION 


Sine series coefficients for defining 
the thin ship 


Wave resistance coefficient of the 2 
model without the thin ship Cw = Rw/(0LV 2 


Wave resistance coefficient of the 2 
model with the thin ship Cww = Rww/(s,L'V>) 


Spline equations defining the non 
optimum thin ship 


Model length, m 


Wave resistance of model without 
the thin ship 


Wave resistance of model with the 
thin ship 


2 
Wetted surface, m 


Model speed, m/s 


Distance along centerline nond imension- 
alized by length (positive aft) ; 


Distance from waveprobe to F.P., 
in meters 


Endpoint of thin ship (nondimesional- 
ized by L) 


Starting poing of thin ship (non- 
dimensionalized by L) 


Distance from model centerline to 
waveprobe along Y axis, in meters. 


Water density 


Thin ship beam at x. 


ENGLISH/SI EQUIVALENTS 


degree (angle) 
foot 

foot per second 
inch 

knot 

lb (force) 

lb (force) - inch 
long ton (2240) 


horsepower 


0.01745 rad (radians) 

0.3048 m (meters) 

0.3048 m/sec (meters per second) 
25.40 mm (millimeters) 


0.5144 m/s (meters per second) 
4.448 N (Newtons) 


0.1130 N-m (Newton-meter) 
1.016 metric tons, or 1016 kilograms 


0.746 kW (kilowatts) 


vi 


ABSTRACT 


The revision to an existing computer program 
which implements wave pattern analysis to improve 
a ship hull form is described. The revision allows 
predictions of the change in wave resistance due to 
adding a given thin ship to a basic hull form. The 
instructions on program usage are included. 


ADMINISTRATIVE INFORMATION 


This Project was authorized and funded by the Naval Material 
Command (NAVMAT) Ship Performance and Hydromechanics Program under 
Program Element 62543N, Subproject Number 43-421-001, Work Unit 
Number 1500-104-32. 


INTRODUCTION 

In 1974 Baba? proposed a method for applying wave pattern analysis to 
the improvement of ship hull forms. His procedure determines a thin ship 
which, when added to a given basic hull form, would reduce the wave making 
resistance. The wave pattern data would be obtained from longitudinal 
wavecut experiments. 

This method was used in an early version of the HULIMP computer 
program. However, the cosine series representation of the thin ship used 
in the original method did not satisfy the prescribed boundary conditions. 
HULIMP was later modified? using a sine series representation of the thin 
ship to satisfy the boundary conditions. 

These early versions of HULIMP would find the optimum thin ship that 
would reduce the wave resistance of a given ship at a given speed and 
thin ship volume. These programs could not calculate the effects on wave 
resistance due to adding an arbitrarily shaped thin ship to an existing 


ship. This omission seriously limited their usefulness as a design tool. 


tReferences are listed on page 5 . 


Because the shape of the optimum thin ship changes with speed, the optimum 


thin ship found at one speed would have to be treated as an arbitrary thin 
ship at another speed. This means that an analysis of the effects upon 
wave resistance due to a given thin ship could not be made through a range 
of speeds. Further, the effects due to design compromises to the optimum 
thin ship cannot be investigated. This new version of HULIMP has the ability 
to find the optimum thin ship for reducing wave resistance and to determine 
the effects on wave resistance due to adding an arbitrary (non-opt imum) 
shaped thin ship. 

Below are additional program documentation covering only the changes to 
HULIMP. The unchanged subroutines are documented in Amato et ait Instructions 


on program usage and a sample output with an arbitrary thin ship are included. 


ARBITRARY THIN SHIP 


The equation of a given thin ship is represented in the program by 


N x-X 
n(x) = L » Cc sin ul) 
n X -X 
n=1 es 


where n(x) is the half beam of the thin ship, N the total number of 
terms in the series, and X, and X, the forward and aft limits of the thin 
ship. The values of c. are derived either from the optimization of the thin 


ship, or from a Fourier sine series analysis of the splines fitted to 


t' 2 half beam offsets of the thin ship. 


When the sine series is fit to a specific thin ship, the equation 


defining the values of c is: 


2 


x, 
a ae i+] x 
A (X,-X) 2 £. (x) sin na aes dx 
i=l i X -X 
at e s > 


where I is the number of spline segments, x; is the longitudinal position 
of the thin ship offset, and f(x) is the cubic spline segment between 
xy and Xe ee These values of Ce are directly substituted into the wave 


resistance equations, replcicing the values of Cc. that correspond to the 


optimum thin ship. 


OPTIMUM THIN SHIP 
The program remains essentially uncuen cd Cor the caiculation of the 


optimum thin ship as described by Amato et si* 


NUMERICAL CONSiDERATIONS 

For the analysis with the optimum thin ship, an anoconriate umber of 
terms is equal to 20 x (XK, - K/L. The exact number will depend upon 
the relative magnitude of the last few terms in the sine series. Excessively 
large numbers of terms may not give better results; it may return values 
of thin ship half-beams in the order of 10000 km. 

For the analysis with a given (non~-optimum) thin ship, up to 30 terms 
can be used. The relative fit of the sine series to the input values are 
printed out. Up to 39 offsets may be input. Even though the cubic splines 
are used for calculating the sine series coefficients, there should be fewer 
terms in the sine series than offsets because the accuracy of the cocflicients 
will be affected by the number of offsets used. The end constraints require 
that the thin ship must have zero beam at each end. 

OVERVIEW OF THE DIGITAL PROGRAM 

Program HULIMP is a FfRTRAN language digital computer program which 
implements Baba's method of ship hull form improvement. Although intended 
to be machine independent, HULIMP has been executed only on the Control Data 


Corporation (CDC) 6000 Series computer installations at DINSRDC. 


A detailed description of the card input data is given in Table 1, and 
a detailed description of the output is given in Table 2. The input data 
consists of a three card input to subroutine READ], a three card input to 
subroutine SUHULL, and available length wave profile input to subroutine 
SUHULL from a card input stream or optional magnetic tape or disk binary 
file. If the non-optimum thin ship analysis is performed, the wave profile 
data are followed by a variable length thin ship beam offset input from 
cards. The last thin ship offset card is followed by a card with 99999. 
starting in column 1. An example input is given in Table 3 and the 
corresponding output in Table 4. 

Run sets can be stacked so that more than one analysis can be performed 


in one job (see MORE, Table 1, card 1) with the option of reading a new 


wave profile or retaining the profile from the previous analysis (see 
REDSUB, Table 1, card 2). 


The wavemaking resistance coefficients presented in the output, Cw 


2 ; 
and Cww, are nondimensionalized using L instead of using the wetted 
2 
surface, S. These values must be multiplied by L/S to convert them to the 
standard wavemaking resistance coefficient form. The beam and draft 


2 , 
constraints have been nondimensionalized by L and L * draft, respectively. 


The job cards required to run the program on the CDC 6000 series 
computer are shown in Table 5. The core requirement for loading and running 
the program is 120,000 octal words. The core requirement can be reduced 
to 57,700 octal words by using segmentation. The reduction in the required 
core size by segmentation significantly increases the execution priority 


of the job. Execution times normally do not exceed 30 seconds. 


Included in the 19 subprograms are 18 operational subroutines and one 
BLOCK DATA subprogram which initializes certain data in particular common 
blocks. The initialization takes place before program execution. 

A block diagram for the digital program is illustrated in Figure l. 

As indicated, the main program, HULIMP, serves as a driver for the entire 
system. In addition to the brief functional description of each subroutine 
given in Table6, a more detailed description of the new subroutines and 
common blocks is provided in Appendix A. 

Standardized communication is used among various subroutines wherever 
possible. FORTRAN symbols used in communicating between routines are defined 
in a glubal sense, i.e., the FORTRAN symbol, DRAFT, is the nondimensional 


draft of the thin ship in all subroutines. 


CDC job control statements, segmentation information, sample innut 


data, and corresponding printer output complete the program description. 


REFERENCES 


1. Baba, E., "Ship Form Improvement by Use of Wave Pattern Analysis," 


Japan Shipbuilding and Marine Engineering, Vol. 8, No. 1, pp 35-43 (1974). 


2. Reed, A.M., and Amato, A.J., "A Ship Hull Form Improvement 


Technique Based on the Baba Theory", DINSRDC/SPD-0820-02, June 1980. 
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TABLE 1 
CARD INPUT DATA DESCRIPTION 


Rec. Format Col. Name Description 
1 14 1-4 NG Number of longitudinal inter- 
vals along hull of existing 
ship 
14 5-8 NK Number of coefficients in 


truncated sine series 


14 9-12 MORE Control variable for process- 
ing more data: 


0 = No more data to be read 


1 = New set of data to be read 
14 13-16 NP1 Control variable for printing 
(if NPl=1) free-wave spectra 
14 17-20 NP2 Control variable for printing 
(if NP2=1) free-wave amplitude 
2 F8.5 1-8 Xs Nondimensional x-coordinate of 
thin ship F.P. 
F8.5 9-16 XE Nondimensional x-coordinate of 
a thin ship A.P. 
F8.5 17-24 DRAFT Draft of thin ship, nondimen- 


sionalized by model length 


F8.5 25-32 REDSUB Control variable to read (if 
REDSUB=0.0) new wave cut 


F8.5 33-40 CS1 (2/L)n(L/2), nondimensional 
beam constraint amidships of 
the thin ship 


F8.5 41-48 cs2 (2/L2) en (x) dx, nondimen- 
x, 
sional volume constraint for 
thin ship 
3 14 1-4 NSC Number of constraints on 
thin ship 


If NSC=3, then the non-optimum thin 
ship case will be analyzed, and 
CS1,CS2, NS1, and NS2 are ignored. 


| 
| 
| 


Table 1 - Card Input Data Description (Cont) 


14 5-8 NS1 Priority of beam constraint 
14 9-12 NS2 Priority of volume constraint 
4 12A6 1-72 TIT Alphanumeric title card 
variable 
5 F8.4 1-8 FN Froude number 
F8.4 9-16 SLWL Length of model, in meters 
F8.4 17-24 Y Transverse location of wave 


cut, in meters 


F8.4 25-32 X Location of first point in 
wave cut relative to origin, 
in meters 


F8.4 33-40 SAMPT Time interval between data 
points, in seconds 
F8.4 41-48 SCL Factor to convert wave input 
to meters 
6 15 1-5 NDATA Number of data points in wave 
profile 
15 6-10 MTR Final truncation point for 
,' the longitudinal wave cut i 
data : 
f 
15 11-15 IOUT Control variable for input of 


wave profile as follows: 


an ee 


0 = Read from magnetic tape 
binary file 


be 
" 


Read formatted data from 
card stream 


15 16-20 MEAN Number of points used to 
compute mean wave height 


I5 21-25 MTRUNC First truncation point fer 
the longitudinal wave cut 
data 


- 
i 
j 
‘ 
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Table 1 - Card Input Data Description (Cont) 


7* 8F10.0 1-80 WAVE(N) Wave height at correspond- 
N=l,.. ing x-y coordinates of 
NDATA longitudinal wave cut - 


(NDATA/8) of these records 
are read (a maximum of 250 
records) 


Bxx F10.5 1-10 XVAL Nondimensional (x/L) distance 
of offset value from the bow. 


F10.5 11-20 YVAL Half beam of the thin ship 
(in meters). The thin ship 
is assumed to be wall sided. 
A maximum of 39 offsets 
(using 39 cards) is allowed. 


Qkx F10.5 1-10 96999, In columns 


1-10 will indicate the end 
of beam offset information. 


* Note: An option exists (see IOUT, record 6) to read wave profile data 
via magnetic tape binary file (tape unit 1). 
ek Omitted if the optimized thin-ship is desired. 


TABLE 2 
OUTPUT INFORMATION 


The line printer output includes: 
o Echo of the first six input records 


o Processed wave cut data (translated vertically to zero mean and 
converted to meters) 


o Necessary information associated with Newman-Sharma method, 
including truncation correction naranetets 


o Free-wave spectra and wave-making resistance of the existing 
ship with and without truncation correction, and the 
influence function of the existing ship 


* o Thin ship half beam offsets echo 
* o Coefficients of the sine series representation cf the input beam values. 
* o Comparison of input half beam values to values calculated “von the sine 


series 


o Coefficients of the sine series (solution vector for the linear 
system) 


o Non-dimensional beam and volume gonstraints 
(Nondimensionalized by L and L¢ * draft, respectively) 


*k o Determinant of the matrix and residuals for the approximate 
solution vector, along with relative and absolute error 


o Wave resistance of the existing ship and modified ship 
(Nondimensionalized using L2, not wetted surface) 
o Influence function for the thin ship and modified shi, 


o Thin ship element station areas and offsets 


o Free wave amplitude before and after modification 
An example of the line printer output is provided in Table 4. 


* Printed only for the non-optimum thin ship case (NSC = 3) 
xk Printed only for the optimum thin ship case (NSC = 3) 


10 


amin 


TABLE 3 


EXAMPLE CASE INPUT 


CARD No. VARIABLE VALUE 
i NG 40 
NK 10 
MORE 0 
NP1 1 
NP2 1 
2 XS -.025 
XE 0.5 
DRAFT 0.0474 
REDSUB 0 
cS1 0 
CS2 0 
3 NSC 3 
NSl 2 
NS2 1 
4 TIT Model 5079 FN=0.28 = 6,521 
FPS RUN 15 
5 FN 0.2793 
SLWL 5.158 
Y 1.364 
m xX 1,075 
SAMPT 0.01 
SCL -,07282 
6 NDATA 1700 
MTR 1679 
IOUT 0 
MEAN 25 
MTRUNC 1552 
8 Thin ship offsets ....... 
9 99999. 


disk (IOUT=0). 


NOTE: Card 7 has been omitted since the wavecut data is stored on the | 
{ 
{ 
| 
{ 
| 


40 10 Q 1 1 


-.025 5 .0474 
ae a 

MODEL 5079 FN = 0.28 C VM = 6.521 FPS RUN 15 

.2793 5.158 1.364 1.075 .01 -.07282 
| 1700 1679 0 25 1552 

-.025 0.0 

0.0 0.0117 

025 60124 

050 0.0124 

.075 0.0109 

-100 .0077 

125 0030 ; 

-150 ~.0018 

2175 ~.0077 

. 200 ~.0131 

-225 ~.0157 

.250 ~.0161 

275 ~.0106 

- 300 ~,0051 

.325 ~.0004 

.350 .001t 

.375 .0018 

.400 0015 

425 0007 

.450 0007 

»475 0 

500 .0 

99999, 


5 
i 
A 
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Table 3 - Example Case Input (Cont) 
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TABLE 5 


Job Control Statements for HULIMP 


The following descirbes the deck arrangement and job control statements 
required to run the program using segmentation. 


JOB CARD 
CHARGE CARD 
ATTACH, HULIMP, ID=XXXX. 
ATTACH, TAPEL ,wAVECUT, LD=CilX™. 
(if the wavecut data are restored in a file) 


SEG:.OAD. 
LOAD(HULIMP) 
EXECUTE. 
7/8/9 END OF RECORD 
HULIMP INCLUDE HULIMP, BLKDAT 
SUHULL INCLUDE SUHULL, TRUNC,WRITE2 ,SUFTRO, THETAX, FRES 
| MATRIX INCLUDE MATRIX,PRIME, COMPUTE 
| READ3 INCLUDE READ3,NEWCOF,SPLVAL,SPINTG,SPFIT 
| READL INCLUDE READ1,DECOMP,SOLVE ,WRITEL 
HULIM TREE HULIUP-(READ1, READ3 ,MATRIX, SUHULL) 
GLOBAL CONST, TRNCWR, COMP ,GLOBL, COMWRZ,HULWRP2, COMWR] 
READ3 GLOBAL RE3COF 
MATRIX GLOBAL MATCOM-SAVE 
SUHULL GLOBAL WAVES 
END 
7/8/9 END OF RECORD 
USER DATA 
: 7/8/9 
6/7/8/9 END OF FILE 


The core requirement for loading and running the program is 57,700 
octal words, using segmentation. If segmentation is not used, 

the core requirement increases to approximately 120,000 octal words 
which significantly lowers the execution priority of the job. Below 
are the job control statements used to run the program without 
using segmentation. 


JOB CARD 

CHARGE CARD 

ATTACH ,HULIMP,, ID=XXXX. 

ATTACH , TAPE1 , WAVECUT , ID=CHXxX. 

HULIMP. 

7/8/9 END OF RECORD 
USER DATA 

7/8/9 END OF RECORD 


6/7/8/9 END OF FILE 


TABLE 6 
FUNCTIONAL DESCRIPTION OF SUBROUTINES 
HULIMP Drives the entire system (calls subroutines), stores 


free wave spectra and computes residuals for solution 
of linear system of equations. 


i 
4 
i 
! 


BLOCK DATA Initializes certain program variables which are in 
common block storage before execution. 

COMPUTE Computes the thin ship offsets and wave-making 
resistance of existing ship and modified ship. 

DECOMP Performs matrix triangularization by Gaussian 
elimination. 

FRES Calculates the value of Fresnel integrals. 

MATRIX Calculates matrix and vector elements for the linear 
system. 

NEWCOF Calls subroutines to calculate the non-optimum thin 


ship sine series. 


PRIME Determines the free wave spectra components due to 
thin ship. 
READ1 Reads and echoes control variables, ship parameters 


and constraint parameters. 
READ3 Reads and echoes the offsets of the thin ship. 


SOLVE Solves a linear system of equations whose coefficient 
matrix has been triangularized. 


SPFIT Fits cubic splines to non-optimum thin ship offsets. 

SPINTG Calculates the integral of f(x)* sin(nx) for the 
non-optimum thin ship. 

SPLVAL Evaluates the splines from SPFIT at a given x. 

SUFTRO Analyses a longitudinal wave cut using the Newman- 


Sharma method, 
SUHULL Reads and checks experimental wave cut data. 


THETAX Computes free wave spectra for the Newman-Sharma method. 
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TABLE 6 (cont inued) 


FUNCTIONAL DESCRIPTION OF SUBROUTINES 


TRUNC Calculates truncation correction parameters by leant 
squares. 

WRITE1 Prints the results of optimization. 

WRITE2 Prints the results of wave pattern analysis. 


— NNER 
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APPENDIX A 


DOCUMENTATION OF ADDITIONAL SUBROUTINES 


PROGRAM HUT IMP 
PURPOSE 


To serve as a driver for the entire program which determines the 


optimum hull torm modification to reduce wave-makiny resistance. 


CALLING SEQUENCE 


None - main routine. 


COMMON BLOCKS 
COMP , COMWR2, CONST, GLOBL, HULWR1,TRNCWR 


SUBPROGRAMS CALLED 
READ1, SUHULL, SUFTRO, PRIME, MATRIX, DECOMP, SULVE, READ3, NEWCOF, 
COMPUTE, WRITE1 


DETAILED DESCRIPTION 
This routine calls subroutines, stores tree wave spectra and computes 


residuals for the solution ot a linear system of equations. 
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SUBROUTINE NEWCOF 
PURPOSE 
To call subroutines that fit a cubic spline to the half beam offsets — 
the non-optimum thin ship and to calculate the respective fourier sine series 
representation. 


CALLING SEQUENCE 
NEWCOF (NTERMS, SLWI.,COFRFS, XS, XE) 


Variable Type Dimension Use Description 

COEFS Real (40) output Sine series coefficient - 
for the non-optimum thir 
ship. 

NTERMS Integer - input Number of terms in the Seri: 

SLWL Real - input Model length 

XS Real - input Nondimensional location 


of the thin ship bow 


XE Real - input Nondimensional location 
of the thin ship stern 


COMMON BLOCKS 


RD3COF 


SUBPROGRAMS CALLED 


SPFIT, SPINTG, SPLVAL 


DETAILED DESCRIPTION 


This subroutine calls subroutines to fit the non-optimum thin ship 
offsets with a cubic spline, and thento find the sine series representat ieu. 


It also prints a comparison of the original thin ship otfsets to the valu.s 


calculated using the sine series. 


SUBROUTINE READ3 
PURPOSE 


To read the half beam offsets of the non-optimum thin ship. 


CALLING SEQUENCE 
READ 3 


COMMON BLOCKS 
RD3COF 
SUBPROGRAMS CALLED 
NONE 
DETAILED DESCRIPTION 
This subroutine reads and echos the half beam offsets of the 


optimum thin ship. 


SUBROUTINE SPFIT 


PURPOSE 
To fit non-parametric spline segments to the non-optimum thin ship 
offsets. 


CALLING SEQUENCE 
SPFIT(X,Y, ELEMS,NPTS) 


Variable Type Dimension Use Description 

ELEMS Real (4,NPTS-1) output Array containing the spline 
information 

NPTS Integer - input Number of points 

X Real (NPTS) input X values 

Y Real (NPTS) input Y values (half beams) 


COMMON BLOCKS 


None 


SUBPROGRAMS CALLED 


None 
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SUBROUTINE SPINTG 


PURPOSE 
To evaluate the integrals f(x)* sin(AX) and f£(x)* cos(AX) where 


f(x) is a cubic spline. 


CALLING SEQUENCE 
SPINTG (XA,XB,X,NPTS, ELEMS,A,CINTG, SINTC) 


Variable Type Dimension Use Description 

A Real - Input Constant for specific 
integral to be evaluated 

ELEMS Real (4,NPTS-1) Input Spline segments generated 
by SPFIT 

NPTS Integer - Input Number of values in 
X-array 

X Real (NPTS) Input X values 

XA Real - Input Lower limit of integration 

XB Real = Input Upper limit of integration 

CINTG Real - Output Integral of £(x)* cos (AX) 

SINTG Real - Output Integral of f(x)* sin (AX) 


COMMON BLOCKS 


None 


SUBPROGRAMS CALLED 
SPLVAL 
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SUBROUTINE SPLVAL 
PURPOSE 


To evaluate a spline at some x value. 


CALLING SEQUENCE 
SPLVAL(X, NPTS, ELEMS,X0,Y0,SO,I ELEM) 


Variable Type Dimension Use Description 

ELEMS Real (4,NPTS-1) Input Spline segments from SPFIT 

NPTS Integer - Input Number of values in x array 

X Real (NPTS) Input X array 

X0 Real - Input X value to evaluate spline 
at 

YO Real - Output f (X0) 

so Real - Output Second derivative at X0 

IELM Integer - Output Index of spline segment 


containing XO 


COMMON BLOCKS 


None 


SUBPROGRAMS CALLED 


None 
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FORTRAN 
Symbol 


XVAL 


YVAL 


NVAL 


Math 
Symbol 


N(x) /2 


COMMON BLOCK RD3COF 


Type 


Real 


Real 


Integer 
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Description 


Nondimensional distance to bow, 
X/L 


Half breadth of thin ship 
corresponding to XVAL 


Number of elements in XVAL and 
YVAL 
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DTNSRDC ISSUES THREE TYPES OF REPORTS 


1. DTNSRDC REPORTS, A FORMAL SERIES, CONTAIN INFORMATION OF PERMANENT TECH- 
NICAL VALUE. THEY CARRY A CONSECUTIVE NUMERICAL IDENTIFICATION REGARDLESS OF 
THEIR CLASSIFICATION OR THE ORIGINATING DEPARTMENT. 


2. DEPARTMENTAL REPORTS, A SEMIFORMAL SERIES, CONTAIN INFORMATION OF A PRELIM- 
INARY, TEMPORARY, OR PROPRIETARY NATURE OR OF LIMITED INTEREST OR SIGNIFICANCE. 
THEY CARRY A DEPARTMENTAL ALPHANUMERICAL IDENTIFICATION. 


3. TECHNICAL MEMORANDA, AN INFORMAL SERIES, CONTAIN TECHNICAL DOCUMENTATION 
OF LIMITED USE AND INTEREST. THEY ARE PRIMARILY WORKING PAPERS INTENDED FOR IN- 
TERNAL USE. THEY CARRY AN IDENTIFYING NUMBER WHICH INDICATES THEIR TYPE AND THE 
NUMERICAL CODE OF THE ORIGINATING DEPARTMENT. ANY DISTRIBUTION OUTSIDE DTNSRDC 
MUST BE APPROVED BY THE HEAO OF THE ORIGINATING DEPARTMENT ON A CASE-BY-CASE 
BASIS. 


